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ABSTRACT 

A tight mass-temperature relation, M{r)/r oc Tx, is expected in most cosmological 
models if clusters of galaxies are homologous and the intracluster gas is in global 
equilibrium with the dark matter. We here calibrate this relation using 8 clusters 
with well-defined global temperatures measured with ASCA and masses inferred from 
weak and strong gravitational lensing. The surface lensing masses are deprojected in 
accordance with N-body simulations and analytic results. The data are well-fit by the 
mass-temperature relation and are consistent with the empirical normalisation found 
by Evrard et al. (1996) using gasdynamic simulations. Thus, there is no discrepancy 
between lensing and X-ray derived masses using this approach. The dispersion around 
the relation is 27 per cent, entirely dominated by observational errors. The next gen- 
eration of X-ray telescopes combined with wide-field HST imaging could provide a 
sensitive test of the normalisation and intrinsic scatter of the relation resulting in a 
powerful and expedient way of measuring masses of clusters of galaxies. In addition, as 
M{r)/r (as derived from lensing) is dependent on the cosmological model at high red- 
shift, the relation represents a new tool for determination of cosmological parameters, 
notably the cosmological constant A. 
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1 INTRODUCTION 

Clusters of galaxies are the largest gravitationally bound 
structures in the Universe and are as such excellent probes 
of cosmic structure formation and evolution. The ensem- 
ble properties of clusters expected in various cosmological 
scenarios can be used to derive constraints on the power 
spectrum of the initial density perturbations and on cosmo- 
logical parameters such as f2o and A (e.g.. Eke, Cole & Frenk 
1996; Bahcall, Fan & Cen 1997; Oukbir & Blanchard 1997; 
Bartelmann et al. 1998; de Theije, van Kampen & Slijk- 
huis 1998). On the scales of individual clusters the inferred 
baryon mass fraction can be used to constrain JIq (White et 
al. 1993; Evrard 1997). In such studies, an important quan- 
tity is the total cluster mass or any observed quantity which 
is tightly related to the mass. 

A promising mass estimator is the mean emission- 
weighted temperature, Tx, of the hot intracluster medium 
(ICM) in clusters of galaxies. Based on numerical simula- 
tions, it has been shown that Tx is a better indicator of the 
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total mass of a cluster than any other optical or X-ray prop- 
erty (Evrard 1990). Recently, Evrard, Metzler & Navarro 
(1996, hereafter EMN) and Eke, Navarro & Frenk (1997) 
showed that there is a tight relation between the mass of a 
cluster and its global X-ray temperature in cosmological gas- 
dynamic simulations, irrespective of the state of the cluster 
(e.g., not restricted to clusters with a 'regular' appearance or 
'isothermal' clusters) and the assumed cosmological model. 
In the simulations, it was found that mass predictions using 
this method (which only involve temperatures) are twice as 
precise as those derived using the /?-model (which require 
the surface brightness distribution in addition, i.e., more 
photons and higher spatial resolution). However, the nor- 
malisation of the relation hinges on numerical simulations 
which may not comprise sufficient detail (EMN; Anninos 
& Norman 1996). Therefore, it is essential to calibrate this 
relation from an observational point of view, by using inde- 
pendent mass estimators. 

The purpose of this Letter is to provide a first observa- 
tional calibration of the M-Tx relation using the relatively 
'clean' way of determining independent cluster masses by 
gravitational lensing. This technique essentially probes the 
projected mass along the line of sight. It is also pointed out 
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that the relation holds the promise of providing a test of the 
geometry of the Universe which is particularly sensitive to 
A. Throughout this Letter, however, we assume a standard 
homogeneous Einstein-de Sitter Universe with Ho = 100ft 
km s"^ Mpc~\ h = 0.5, flo ^ 1 and A = 0. 



2 THE MASS-TEMPERATURE RELATION 

Navarro, Frenk & White (1997, hereafter NFW) found in 
their numerical simulations that the dark-matter distribu- 
tion in present-day clusters have self-similar density profiles 
when the radial coordinate is scaled to the radius contain- 
ing an overdensity of 5 = 200 relative to the critical density. 
More precisely, defining the overdensity as 



S{rs, z) 



47r/9c(z)r-f 



(1) 



where pc{z) = Pco(l + zf and pco = 3Jfo/(8vrG), NFW 
found that clusters are well-described by the density profile 
p{x) oc x~^{l -I- cx)~'^ , where x = r/r2oo, in any cosmol- 
ogy. The variation in c « 5-10 with mass, cosmological pa- 
rameters, and redshift is small (Cole & Lacey 1996; NFW; 
Bartelmann et al. 1998; Eke, Navarro & Frenk 1997) and so 
clusters form a homologous family to a good approximation 
when scaled to a given overdensity. Optical (Carlberg et al. 
1997) and lensing observations (Fischer & Tyson 1997) seem 
to support this conclusion. 

For a cluster in quasi-equilibrium (Natarajan, Hjorth & 
van Kampen 1997) the virial theorem for the dark matter 
states that M{r) cx r (v^^ . Self-similarity implies that the 
constant of proportionality depends on the adopted overden- 
sity only. Finally, the global X-ray temperature is assumed 
to be proportional to the global mean velocity dispersion 
of the dark matter (at any time), i.e., Tx oc (v'^) • For 
example, this would be the case in the absence of tran- 
sient effects and non-gravitational heating or cooling ef- 
fects. In the case of equipartition one would have a uni- 
versal f3 = fiTUp {y^^ /{kT) — 1. Combining these assump- 
tions (quasi-equilibrium, self-similarity, proportionality be- 
tween dark-matter velocity dispersion and X-ray tempera- 
ture) leads to a simple scaling relation between the charac- 
teristic mass and radius at a given overdensity, S, and the 
global emission-weighted temperature of the hot X-ray gas. 



Ms{rs) = ksrgTx, 



(2) 



where kg is a constant depending on 5. This equation ex- 
presses the structural invariance of clusters under mass and 
redshift transformations and does not rely on any particular 
dark-matter density profile or the assumption of hydrostatic 
equilibrium. 

Combined with the definition of the overdensity (eq. 1) 
this expression leads to the mass-temperature relation, 



Ms 



t.3/2 
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or, equivalently, the size-temperature relation, 
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Rather than trying to compute the prefactor kg from 
first principles, EMN used numerical simulations to calibrate 
these relations. They found the radius rsoo to be a conserva- 
tive estimate of the boundary between the virialised region 
of the clusters and their outer envelopes. At z — 0.04 using 
5 = 500 they found a universal prefactor independent of flo, 



Msoo = 2.22 X 10^^ 



Tx 



10 keV 



3/2 



Mr. 
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In the simulations the scatter around this relation was found 
to be only 15 per cent compared to 30 per cent when using 
the /3-model to estimate the mass. 

Mohr & Evrard (1997) have recently shown that obser- 
vations of nearby clusters lead to an intrinsic scatter of 10-15 
per cent in the relation between cluster isophotal size and 
mean emission-weighted temperature Tx (similar to eq. 4) 
regardless of the state of the cluster (merging, cooling fiow) 
thus giving added support to the existence of a tight mass- 
temperature relation. 



3 LENSING MASSES 

In order to test and calibrate the mass-temperature relation 
observationally we shall use independent masses determined 
from gravitational lensing. Since lensing masses are given in 
the literature as a function of physical radius rather than 
overdensity, we shall use eq. (2) to express the temperature 
as a function of M{r)/r instead of Ms - Thus the relation we 
shall test observationally is 



M{r) 



1015 M, 
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1 Mpc 



— kg 



Vio kevy 
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where rs{l+z) is the co-moving angular radius of the cluster. 



For an isothermal sphere, M{r) oc r, kg would be a constant 
independent of radius or overdensity. However, given the fact 
that clusters are described by more complicated density and 
temperature profiles (NFW; EMN) ks varies slightly with S 
in the range considered. While an overdensity of ~ 500 was 
recommended (rsoo ~ 1-2 Mpc), EMN provided normali- 
sations for S = 100, 250, 500, 1000, 2500. Converting these 
into equivalent values for kg we find 0.76,0.91,1.01,1.09,1.14. 
These slowly varying numbers are used to compute kg as a 
function of 5 by spline interpolation. 



3.1 Deprojection 

Lensing provides the 2D projected (surface) mass, M2u{R) 
(where R indicates a projected radius), of the cluster. In 
general the 2D mass at a given radius is larger than the 
3D mass, M3o{r), evaluated at the same radius (r = R). 
The best way of obtaining a deprojection relation for M30 
(which is the quantity entering eq. 6) is to study numerical 
simulations of galaxy clusters, preferably a fair sample of 
these. We have used the catalogue of simulated standard 
CDM clusters (fi = 1) of van Kampen & Katgert (1997) 
to find such a relation (van Kampen, in preparation). As 
the clusters we shall study in this Letter are biased towards 
massive clusters, we selected only clusters with a total mass 
within the Abell radius (3 Mpc) of at least 10^^ M© . The 
deprojection relation M3u{x)/M2d{x) for these 41 clusters 
is plotted as a function of dimensionless radius x = r/r2oo in 
Fig. 1, along with its scatter, which is fairly substantial at 
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Figure 1. This figure shows M3d(3;)/M2d{2;) as a function of 
X = r/r200 for the Hernquist (1990) profile (dashed curve) and the 
model of BBS (dashed-dotted curve). The solid curve is the corre- 
sponding mean deprojection factor for simulated massive clusters 
in a CDM O = 1 Universe (see text for details) and the dot- 
ted curves indicate the la confidence interval. The corresponding 
curve for an open CDM = 0.2 Universe (not plotted here) lies 
10 per cent higher in good agreement with the Hernquist model. 

small radii as substructure along the line of sight becomes 
important for the projected mass. In an open fJo = 0.2 CDM 
Universe the corresponding curve (not plotted) is 10 per cent 
higher due to the smaller influence of substructure. 

In Figure 1 we also plot the deprojection relation of 
the Hernquist (1990) and the Brainerd, Blandford & Smail 
(1996, hereafter BBS) models which both have analytic de- 
projection properties (the BBS model is the limit of 77 ^ oo 
of the Hjorth & Kneib (1998) model). As the total mass of 
the NFW model is infinite it is not useful for this purpose. 
However, we have made use of the fact that the half-mass 
radius in the Hernquist model is roughly equivalent to r2oo 
in the NFW model for c « 5 (Cole & Lacey 1996). In the 
outer parts the Hernquist and BBS models coincide, but 
in the centre there is a marked difference between the two 
curves because of their differing divergence properties. The 
Hernquist model which diverges as p ~ in the centre 
is similar to the NFW profile and M3D(a:)/M2D(a:;) for 
X — > while the BBS model, which has a stronger central 
cusp p ~ r~^, tends to the value for the singular isother- 
mal sphere 2/7r ~ 0.64. This shows that deprojection of 2D 
masses at small radii depends sensitively on the exact slope 
of the inner cusp of dark-matter density profiles (Fukushige 
& Makino 1997; Moore et al. 1998; Kravtsov et al. 1998). We 
finally note that the Hernquist model is in excellent agree- 
ment with the numerical results of the open model. 

In this Letter we shall use the relation as a function of 
proper radius R to deproject the lensing masses. For this 
purpose we introduce a convenient fitting function, 

^{R) = 0.56 tan-^ ( — , (7) 
hhu \^0.28Mpc y ' ^ ' 

where the coefficients have been determined from a non- 
linear least-squares fit up to i? = 2 Mpc. We note that if 



such a deprojection correction is not applied, lensing (2D) 
masses will be higher than X-ray (3D) masses by a factor of 
1.5 on average. 



4 DATA 

We have compiled a list of clusters with well-determined 
X-ray temperatures and masses determined independently 
using gravitational lensing. The data are shown in Table 1. 

The X-ray data used here are from a recent compila- 
tion of temperatures of intermediate and high redshift clus- 
ters observed by the Advanced Satellite for Cosmology and 
Astrophysics (ASCA) (Mushotzky & Scharf 1997). The tem- 
peratures were measured in a uniform way out to a radius 
of 3-6 arcmin depending on the redshift of the cluster. 

The lensing masses are from various studies of individ- 
ual clusters, mostly using the 'weak lensing' method pio- 
neered by Kaiser & Squires (1993), but also from variations 
in number counts of background galaxies (Broadhurst, Tay- 
lor & Peacock 1995; van Kampen 1998). We included only 
clusters with masses determined out to radii larger than 
400 kpc to minimise deprojection and substructure effects 
from the central regions of the clusters. One cluster mass 
(MS 1358+62) was determined from wide-field HST imag- 
ing. The normalisation of the mass of A2163 was adjusted in 
comparison with X-ray masses (derived from the (5 model), 
i.e., this mass is not completely independent of the temper- 
ature (Squires et al. 1997). 

We show the results for the eight clusters in Fig. 2 in 
which we plot M^oiR) / R derived from lensing studies as a 
function of ksTx- The relation predicted by eq. (6) gives an 
excellent fit to the data. The best-fit line has a normalisa- 
tion which is 88 per cent of that predicted by EMN and the 
dispersion (rms) about the relation is 27 per cent in mass, 
somewhat smaller than that expected from the quoted ob- 
servational errors alone. 



5 DISCUSSION 

The observed scatter around the mass-temperature relation 
(eq. 6) is dominated by observational errors and is consistent 
with having no intrinsic scatter. The data thus support the 
existence of an M-Tx relation as a fairly accurate indepen- 
dent estimator of cluster masses. It is, however, important 
to point out that the data set presented here may be affected 
by systematic errors in both the masses and the tempera- 
tures. Adopting different data sets could lead to significant 
changes in e.g. the normalisation of the mass-temperature 
relation. 

X-ray temperatures of hot clusters are usually uncer- 
tain due to the few photons detected above 8 keV with 
ASCA and quoted errors normally do not incorporate possi- 
ble systematic errors. Thus, many compilations (e.g., Sadat, 
Blanchard & Oukbir 1998) may be affected by the fact that 
temperatures often differ from author to author (Mushotzky 

6 Scharf 1997; Allen 1998; Yamashita 1997) due to differ- 
ences in data analysis and use of data from other satellites 
{ROSAT, Gmga). For example, Allen (1998) has shown that 
cooling fiows may bias global temperatures as derived by 
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Table 1. Observational data on clusters with ASCA temperatures and lensing masses. The temperatures are from 
Mushotzky & Scharf (1997). The lensing masses are generally taken from the most recent publication of a given cluster. 
For the data of Smail et al. (1995) we have assumed an uncertainty of 40 % in the masses and adopted the no-evolution 
model for the redshift distribution of faint background galaxies. For the data of Squires et al. (1996ab,1997) we have 
estimated the masses inside 210". All masses and radii are computed assuming h = 0.5, Qq = 1 and A = 0. 



Cluster 


z 




Tx 
(keV) 


R 

(Mpc) 


AhoiR) 

(1014 Mq) 


M-ioiR) 
(1014 Mq) 
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ks 


Reference 


Abell 2218 
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0.80 


9.4+1.7 


6.5+1.2 


2714 


1.14 


Squires et al. (1996a) 


Abell 1689 





18 




0.48 


10.0+1.8 


5.8+1.1 


11057 


1.15 


Taylor et al. (1998) 


Abell 2163 
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19 7+2.0 
-^^■'-2.0 
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13.0+10 
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1.14 


Squires et al. (1997) 


Abell 2390 
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8.901°:?? 
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10+4 
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1.14 


Squires et al. (1996b) 


MS 1455.0+2232 
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r 45+0.29 


0.45 
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3859 


1.14 


Smail et al. (1995) 


MS 1358.4+6245 





33 


6 ^O+O *^* 
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1.00 
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3.2+0.4 
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1.00 


Hoekstra et al. (1998) 


RX J1347-1145 





45 


11.371^-0 


2.00 


34+8 


27+6 


385 


0.98 


Fischer & Tyson (1997) 


MS 0015.9+1609 





54 


8.0tl:0 


0.60 


8.5+3.4 


5.4+2.2 


2355 


1.14 


Smail et al. (1995) 




kJ,(keV) 

Figure 2. This figure shows the M^aiR) / R-Tx relation for lens- 
ing clusters of galaxies. Filled circles are ground-based data, the 
open diamond is the HST data point of MS 1358+62, the open tri- 
angle is A2163. The error bars do not include deprojection errors. 
The solid line is the relation predicted by eq. (6) as normalised 
by EMN and has not been fitted to the data. The dashed line 
minimises the mean relative residual and has a normalisation 12 
per cent lower. 



Mushotzky & Scharf (1997) downward by about 30 % on av- 
erage. In this Letter we have used the Mushotzky & Scharf 
(1997) data because of the uniformity and simplicity of the 
analysis and the possibility of extending it to fainter and 
higher redshift clusters. 

Ground-based lensing data are affected by fairly large 
seeing corrections. HST results are therefore preferable but 
usually only available inside a small radius. The normali- 
sation of weak-lensing masses (due to the mass-sheet de- 
generacy) can be carried in various ways, using e.g. mul- 
tiple arcs, magnification bias, fits of analytic models or 
comparison with X-ray profiles. Sadat et al. (1998) found 
that a normalisation of ks about 64 per cent of the EMN 
value (no deprojection was applied) was consistent with the 
temperatures and masses of 5 HST clusters studied out to 
400 kpc by Smail et al. (1997). Such a small normalisation 
would imply that the data presented here have systemati- 



cally overestimated masses or underestimated temperatures. 
The high-quality HST data point for MS 1358+62 may in- 
dicate that ground-based masses are indeed overestimated. 
However, the mass of MS 1358+62 as derived from weak 
lensing could also be underestimated e.g. due to deprojection 
errors arising because of the high ellipticity of the cluster. 
Such an underestimate of the 3D mass is supported by the 
disagreement between the velocity dispersion derived from 
weak lensing (780 ± 50 km s~^) and that found from direct 
spectroscopic measurements as well as from strong lensing 
(~ 1000 km s"^) (Hoekstra et al. 1998). //ST clusters may 
also have underestimated masses due to the fact that a fit 
of an isothermal sphere to the mean tangential shear inside 
a small radius (Smail et al. 1997) in general biases masses 
low (van Kampen & Hjorth, in preparation). 

It is a long standing discussion whether masses deter- 
mined from lensing agree with or exceed X-ray masses de- 
termined using the /? model (see e.g. Smail et al. 1997; Allen 
1998). If we take the results presented in Fig. 2 at face value 
the good agreement between the EMN normalisation and 
the observational calibration indicates that there is no such 
discrepancy when using the mass-temperature relation. If 
anything the X-ray masses computed using the EMN nor- 
malisation are slightly higher (by ~ 10-20 per cent) (cf. Fig. 2 
and MS 1358+62) than lensing masses. Such an effect would 
be consistent with the predictions of simulations incorporat- 
ing the effects of galactic winds (Metzler & Evrard 1998) 
which contribute to heating the ICM. 

Besides its use as a straightforward mass-estimator for 
any cluster with a well-determined global temperature the 
mass-temperature relation holds the promise of becoming 
an important cosmological tool, bearing a resemblance with 
the Fundamental Plane or TuUy-Fisher scaling relation for 
elliptical or spiral galaxies, respectively, in that it relies on 
simple scaling relations with 10-20 per cent scatter, but pre- 
sumably involves much smaller evolutionary corrections. A 
direct cosmological application of the M-Tx relation would 
be to examine the inferred deviations from it as a function of 
redshift. A possible trend with redshift could be indicative 
of (i) evolutionary effects (ii) the assumed redshift distri- 
bution of the faint background galaxies, N{z), or (iii) the 
parameters entering the assumed cosmological model. 

Typical evolutionary effects could be non-gravitational 
heating or cooling of the ICM such as effects of feedback 
mechanisms like galactic winds which introduce systematic 
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structural changes of the ICM (Metzler & Evrard 1998) or 
coohng flows (AUen 1998). Possible 'outliers' from the rela- 
tion could be due to e.g. merging clusters with a very un- 
settled temperature distribution (Schindler 1996) or highly 
elongated clusters which give large deprojection uncertain- 
ties depending on the viewing angle. 

While the inferred lensing masses of low and interme- 
diate redshift clusters are fairly insensitive to the assumed 
median redshift of the background galaxies, high-redshift 
clusters are very sensitive to the assumed median redshift 
(Small et al. 1995; Luppino & Kaiser 1997) and so devia- 
tions from the expected relation at high redshift could be 
used to constrain N{z). 

Finally, the world model enters through the derived 
masses and sizes via the expression for the angular diam- 
eter distance. In the simplest form, Tx oc Ds/Dls, where 
Tx is a directly measurable intrinsic quantity and Ds/Dls 
is the ratio between the source and lens-source angular di- 
ameter distances. Thus, the method can be used as a test 
for the geometry of the Universe, which is less sensitive to 
inhomogeneities along the line of sight than small standard 
rods/candles (e.g., SN la) (Hadrovic & Binney 1997). In- 
dividual massive high-redshift clusters could therefore be 
fairly unbiased discriminators between different cosmolog- 
ical models, particularly sensitive to the cosmological con- 
stant A. At 2 = 1 the difference between a (f2o, ^a) = (1, 0) 
and a (iloif^A) = (0.2,0.8) Universe is 25 per cent in 
Ds/Dls- Moreover, if the measurement of the global X-ray 
temperature is supplemented with spatially resolved X-ray 
imagery additional constraints on fio can be derived (Sasaki 
1996; Pen 1997). 



6 CONCLUSION 

Based on numerical simulations (EMN; Eke et al. 1997) and 
observations of nearby clusters (Mohr & Evrard 1997) the 
existence of a tight mass-temperature relation has been sug- 
gested. The results presented here provide support for this 
assertion and indicate that the mass-temperature relation 
(eq. 6) can be used to determine cluster masses with a pre- 
cision of 27 per cent (Fig. 2). There seems to be no significant 
discrepancy between deprojected lensing masses and masses 
derived from X-ray temperatures, using the normalisation 
found in numerical simulations (EMN). 

The origin of this tight relation is believed to be the 
fairly simple physics entering the relation (cf. Sec. 2), namely 
virialisation of gravitationally bound structures with self- 
similar dark-matter density distributions that are in global 
quasi-equilibrium with the hot ICM, independent of the cho- 
sen world model, power spectrum, or exact formation red- 
shift of the cluster. 

We have cautioned that the observational data dis- 
cussed in this Letter are quite uncertain and possibly af- 
fected by systematic errors. The results should therefore only 
be taken as an indication of a tight mass-temperature rela- 
tion. However, the future observational situation is promis- 
ing. A sample of clusters with very precise lensing masses 
(e.g., from wide-field HST imaging with the ACS) to about 
10 per cent or better (e.g., Natarajan et al. 1998; Hoekstra 
et al. 1998) and equally accurate temperatures (e.g., with 
AXAF, Spectrum-XG, or XMM) would allow us to study 



the intrinsic scatter of the relation and determine a precise 
normalisation. This could provide a direct and reliable mass 
estimator for distant clusters with important cosmological 
implications. 
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